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The most sensitive experiment to date to search for the muon and electron lepton number violating decay 
K1 — > /i^e^ has detected no events consistent with this process. Based on this result, the 90% confidence level 
upper limit on the branching fraction is B[K% /i^e^) < 4.7 x 10~^^. 
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This Letter reports a significantly improved ex- 
perimental upper limit on the rate for the decay 
K1 ^^e^. This process would violate the 
conservation of muon and electron lepton num- 
ber (referred to as separate lepton number) while 
conserving total lepton number; it is not allowed 
in the Standard Model of particle physics. In- 
corporating neutrino masses and mixing into the 
Standard Model, consistent with current informa- 
tion on these quantities, leads to predicted rates 
for — > fi^e^ well below experimental sensitiv- 
ities Q. Hence, observation of this decay would 
signal new physics processes. 

Experiments to search directly for separate lep- 
ton number violation have been performed for 
many years, all with null results. Some of the 
best limits come from previous searches for Kl — > 
/i^e^ l^jH, for which the combined upper limit 
is 2.4 X 10^^^, and from searches for — > 
^V^eT |, K+ ^ TT+^i+e- |, M+ ^ e+j §, 
/i+ ^ e+e+e- §, and e-N g. The sen- 

sitivity of these processes to mechanisms which 
allow for separate lepton number violation varies. 



Several theoretical models allow 



eluding some at rates as large as current exper- 
imental limits: horizontal gauge interactions 



left-right symmetry ||l|,|lO|, technicol or |ll[ , com- 
positeness |lj], and supersymmetry 13|. 

The experiment (E871) was performed in 
the B5 beam line of the Alternating Gradi- 
ent Synchrotron at Brookhaven National Labora- 
tory (BNL). Data were collected during running 
periods in 1995 and 1996. A 24 GeV proton beam 
incident on a 1.4 interaction length platinum tar- 
get produced the neutral beam. A targeting an- 
gle of 3.75° was chosen to maximize the yield 
while minimizing the n/ ratio. Sweeping mag- 
nets downstream of the target removed charged 
particles. Thin lead foils in the first magnet aper- 
ture reduced the 7 flux in the beam. Three col- 
limators defined the neutral beam to subtend a 
solid angle of approximately 4x16 mrad^. The 
proton intensity was typically 1.5 x 10^^ in a 1.2- 
1.6 s pulse every 3.2-3.6 s, resulting in 2 x 10 
per pulse {2 < pK < 16 GeV/c), 7.5% of which 
decayed in an evacuated decay volume between 
9.75 and 20.75 m from the target. We measured 
the n/K^ ratio to be 8 ± 3. 

The experimental apparatus, illustrated in 
Fig. |l|, has been described previously [|l^. Ac- 
cordingly, we describe here only the most impor- 
tant features, emphasizing characteristics of the 
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detector most germane to this search. A magnetic 
spectrometer consisting of six pairs of tracking 
chambers and two magnets was used to measure 
kinematic parameters of charged decay prod- 
ucts. A beam stop was placed in the first 
magnet to absorb the neutral beam. The spec- 
trometer was followed by scintillation hodoscopes 
(TSC) at two longitudinal positions, used to se- 
lect events with two charged particles and de- 
fine the event time. Both hodoscopes contained 
one X (horizontal) measuring plane on either side 
of the detector. The downstream module con- 
tained an additional plane on either side provid- 
ing measurements in the y (vertical) view. A seg- 
mented threshold Cerenkov detector (CER) and 
lead glass calorimeter (PBG) were used to iden- 
tify electrons. A 30.5 cm thick iron filter fol- 
lowed the PBG. Muons were identified by scin- 
tillation hodoscopes (MHO) and a range finder 
(MRG) downstream of it. Six MHO planes at 
five different z locations provided position and 
timing information. The MRG consisted of 26 
vertical and 26 horizontal planes of proportional 
tube hodoscopes, located between absorbers of 
steel, marble, and aluminum. The amount of ma- 
terial between successive planes increased with 
depth and corresponded to 5% increments of 
muon range. 

The lowest level trigger (LO) required a coinci- 
dence of signals in all six TSC planes, providing 
two X and one y measurements on each side of the 
detector. Further, a parallelism condition was im- 



posed: each pair of x view signals was required 
to be consistent with a particle trajectory with 
\dx/dz\ < 31 mrad. This requirement maintained 
good acceptance for two-body decays while reduc- 
ing the acceptance for the dominant three-body 
modes. The typical LO trigger rate was approxi- 
mately 70 kHz. 

The next trigger level (LI) was designed to ac- 
cept all dilepton decay modes; it was formed from 
a coincidence of an LO trigger and signals from 
particle identification (PID) detectors. The muon 
signal was taken from the MHO plane located at 
a z position corresponding to an energy loss of 
1 GeV. The electron signal was provided by the 
Cerenkov detector. Spatial correlation between 
TSC and PID signals was required. In addition 
to the dilepton modes, one of every 1000 LO trig- 
ger events was selected to form the "minimum 
bias" sample. It was used for detector calibration 
and for flux determination, based on the number 
of K1 — > 7r+7r~ events in that sample. The LI 
trigger rate averaged about 7 kHz. 

Events satisfying the LI trigger were digitized 
and transferred to a set of eight processors in 
which a software trigger (L3) was implemented. 
It did fast event reconstruction using hits in the 
TSC and all tracking detectors. At least one track 
on each side of the spectrometer and a decay ver- 
tex within the neutral beam were required. In ad- 
dition, a two-body invariant mass M+_ exceeding 
460 MeV/c^ and two-body transverse momentum 
Pt less than 60 MeV/c were required for the iie 
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trigger. Charged particle mass assignments were 
determined by the triggering PID detector. An 
event with more than one LI dilcpton hypothesis 
was accepted if any of the triggered modes satis- 
fied the L3 criteria. The L3 trigger algorithm was 
rnn on minimum bias events but no selection on 
kinematic quantities was made. All events pass- 
ing the L3 trigger were written to tape for off-line 
processing. 

The off-line pattern recognition software used 
a similar but more thorough algorithm for pat- 
tern recognition. More rigorous selection criteria 
were applied to geometric and kinematic quan- 
tities. Events were selected if M^- exceeded 
470 MeV/c^ and if either of the requirements 
Pt < 40 MeV/c or 6c < 4.5 mrad was satisfied, 
where 6c is the angle between the kaon direction 
(determined from the target and vertex position) 
and the direction of the momentum sum vector 
of the two charged particles. 

The candidate sample was se- 

lected from events satisfying these criteria and 
having a LI hypothesis. Projected track po- 
sitions were required to be spatially consistent 
with signals in the PID detectors. A particle 
was identified as an electron if it had appropri- 
ate signals in the CER and PBG, the track time 
determined from the tracking spectrometer was 
consistent with the CER time, and the momen- 
tum was consistent with the energy deposited in 
the PBG. A particle was identified as a muon if it 
projected to hits in the MHO and MRG, the track 
time was consistent with the MHO signal times, 
and the energy inferred from the MRG exceeded 
80% of the momentum. Also, any PBG signal as- 
sociated with a muon candidate was required to 
be inconsistent with that of an electron. 

The best estimates of kinematic parameters of 
all events satisfying the above criteria were deter- 
mined using two independent algorithms that had 
different sensitivities to errors in pattern recogni- 
tion. One algorithm (FT) minimized a for 
a kinematic fit to all hits in the tracking detec- 
tors, appropriately accounting for detector reso- 
lution and multiple scattering using an error ma- 
trix. The second algorithm (QT) did separate 
calculations of kinematic quantities for the front 
and back halves of the spectrometer and used the 



two measurements to form a single momentum 
and a measure of the track quality. The FT algo- 
rithm had a better mass resolution (L13 versus 
1.26 MeV/c^ for QT for if£ ^ tt+tt" events) 
and was used to determine the event kinematics. 
It was required that M^e and pt determined by 
the fitters be consistent. 

Events were required to be consistent with orig- 
inating from a single kaon decay. The times of 
the two particles as measured by the spectrom- 
eter were required to be consistent within mea- 
surement uncertainty and their trajectories were 
required to project to a common vertex with 
2; > 9.75 m and transverse coordinates within the 
neutral beam. 

Events were selected on the basis of the qual- 
ity of the kinematic fit in order to reject those 
with tracking mismeasurement or with pion de- 
cay or large angle scattering within the spectrom- 
eter. In particular, the measurements of front 
and back momenta were required to be consistent, 
and events with poor tracking were rejected. 
High momentum particles were particularly sus- 
ceptible to errors in track reconstruction, and a 
maximum momentum of 8 GeV/c was imposed. 
This requirement also ensured that muons were 
contained in the range-finder and that Cerenkov 
signals were not caused by pions. Events were 
rejected if particle trajectories projected to any 
material in the spectrometer volume other than 
the relevant detectors. 

To ensure that event selection criteria were 
free of bias from knowledge of potential signal 
events, they were chosen by studying jjL^e^ 
candidates with M^e > 485 MeV/c^ and p^ < 
900 (MeV/c)^, but excluding a region, 490 < 
M^e < 505 MeV/c^ and p^ < 100 (MeV/c)^, 
larger than the potential n^e^ signal re- 

gion. This exclusion region corresponded to a 
±5 a interval in mass and a 3 cr interval in p^. 

The primary source of — > /i^e^ background 
is K'^ — > ireiy in which a pion decays upstream 
of the muon filter (about 4% of all Kl decays). 
Misidentification of the pion as a muon results 
in a maximum value for Af^e of 489.3 MeV/c^ if 
track momenta and directions are correctly mea- 
sured and the pion is assigned a muon mass. 
The resolution in M„e inferred from the measured 
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Figure 2. Data and Monte Carlo distributions 
of (a) momentum asymmetry and (b) p^}p. Events 
shown satisfy all selection criteria except those on 
the quantities displayed above and have 
493 MeV/c2 and 100 < < 900 (MeV/c)^ 
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K1 TT+TT^ mass resolution was 1.38 MeV/c^; 
hence, background arising from Gaussian tails 
in the mass measurement is negligible, but non- 
Gaussian effects could be important. One source 
of non-Gaussian errors is elastic scattering in the 
vacuum window or in the first tracking detec- 
tor (0.12% and 0.23% radiation lengths, respec- 
tively), referred to as upstream scatters. Scatter- 
ing and decay in the plane normal to Pe x (the 
Kl decay plane) may increase the ire opening an- 
gle and hence the value of M^e, while scattering 
or decay out of the plane cannot increase M^e 
significantly. Monte Carlo simulations showed 
that the dominant background occurs when a low 
energy electron elastically scatters upstream and 
the pion decays upstream of the spectrometer. 

To study this background in detail we imposed 
the requirement Pe > 1 GeV/c and examined the 
remaining high mass events outside of the exclu- 
sion region. The momentum asymmetry and the 



component of the fie transverse momentum in the 
decay plane (denoted by p^) are shown in Fig. ^ 
together with the Monte Carlo predictions. The 
sign of is taken to be positive if it lies on the 
electron side of the spectrometer. The upstream 
scatter events are characterized by large momen- 
tum asymmetry and more importantly by large 
p^}p. Without additional selection criteria, our ex- 
pected background in an appropriate signal re- 
gion would be about one event. We imposed the 
requirement that {p^ —pe)/{Pfi + Pe) < 0.5 and 
that p^ be small ||l^, which reduced this back- 
ground significantly. 

A second potential source of background is ac- 



cidental coincidences of 



nev and 



TTlxv decays. Because the muon and electron orig- 
inate from independent decays, they can recon- 
struct with a value of M^e > Mk- Monte Carlo 
simulations were done to study this background. 
In many of these events at least one of the pion 
trajectories is fully contained in the spectrome- 
ter, and the background results when the x and 
y view tracks on the pion side are mispaired. 
Rejecting events with three or more fully recon- 
structed tracks in the spectrometer reduced this 
background by an order of magnitude. 

The fi^e^ selection criteria, including 

choice of the signal region, were determined by 
simultaneously varying values of relevant event 
selection parameters to maximize the sensitivity 
to signal while suppressing the expected contri- 
bution from the dominant source of background 
to 0.1 event. The calculated background after ap- 
plication of all selection criteria except those on 
M^e and is compared to the data in Fig. ^. 
For M^e < 490 MeV/c^, the background is dom- 
inated by correctly measured — > neiy decays. 



For Mf,e > 



493 MeV/c^, events passing the se- 
lection criteria are dominated by large upstream 
scatters, with a lesser contribution from acciden- 
tals. Reducing the background by an additional 
factor of 10 would require tighter selection criteria 
which would result in a 50% acceptance loss. 

For A/^e < Mk, the signal region was defined 
by the ellipse [pf,/(20 (MeV/c)2)]2 + [AM/(2.4 
MeV/c2)]2 < 1, where AM = M^e - Mk- For 



M^e > Mk, the signal region was defined by 
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Figure 3. Data and Monte Carlo distributions 
of M^e for events with p| < 20 (MeV/c)^ (the 
p"^ range of the signal region). We estimate the 
uncertainty in the Monte Carlo normalization to 
be 10%. The calculated ^±e=F signal curve 

assumes a branching fraction of 2.1 x 10^^^. 
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Figure 4. Plot of versus M^e- The exclusion 
region for the blind analysis is indicated by the 
box. The signal region is indicated by the smaller 
contour. 



AM < 4 MeV/c2 and < 20 (MeV/c)^. The 
background is larger in the region M^e < Mk; 
hence different shapes were chosen for the signal 
region below and above Mk as a compromise be- 
tween acceptance and background rejection. Af- 
ter all selection criteria (including the choice of 
the signal region) were determined, all data (in- 
cluding those in the exclusion region) were rean- 
alyzed. Figure ^ shows the final distribution in 
p^ versus M^e- There are no events in the signal 



fjL^e^ sensitivity is determined from 



region. 

the number of tt+tt^ decays in the mini- 

mum bias sample. These events were required to 
satisfy an appropriate subset of the final selec- 
tion criteria discussed above, and were required 
to have no PBG signals consistent with those of 
an electron. A fit in the p^ versus A/tttt plane was 
done to subtract residual — > ir^v background 
and to determine the number of K'^ —> tt+tt" 
events. Small differences in geometric acceptance 
and cut efficiencies were determined by Monte 
Carlo simulation. 

In the absence of any signal events, the 90% 
confidence level upper limit on the — > fi^e^ 
branching fraction is given by 



2.3 ^B-j^-j^ 



/tttt ^tttt 1 1 6 



PID 



M£ fie ^fie ^fie 



where B^^ is the —>■ TT+TT branching 
fraction, R is the tttt prescale (a hardware factor 
of 1000 times a software factor of 20), Nt^-^ is the 
number of tt+tt" events in the prescaled 

minimum bias sample (including a 0.05% cor- 
rection for Kl ~ Ks interference), /,r7r is a fac- 



tor to account for loss of 



events 



due to hadronic interactions in the spectrome- 
ter, Atjt, and A^e are the mode dependent accep- 
tances (including geometric acceptance and selec- 
tion criteria efficiency), ej^g and are the effi- 



ciencies of the LI and L3 triggers, and e^^^ and 
ej^™ are the efficiencies of the particle identifica- 
tion. The geometric acceptance for — > ^'^e^ 
{Kl vr+vr") decays with 9.75 < z < 20.75 m 
and kaon momentum 2 < pK < 16 GeV/c was 
2.36% (2.63%). The parallelism requirement re- 
duced this to 1.97% (2.21%). Event selection 
criteria further reduced the acceptance to 1.14% 
(1.62%). Table| summarizes the factors entering 
into the B{Kl ^ n^e^) upper limit calculation. 

The resulting 90% confidence level upper limit 
on the branching fraction is B{K^ M^e^) < 
4.7 X 10^^^. This represents the most sensitive 
/i^e^ to date. 



search for A'" 
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0.002067 ±0.000035 




1.14% ±0.006% 


/tttt 


0.959 ±0.0058 


LI 

^ae 


0.974 ±0.0046 


R 


2 X IQ-* 


L3 
^iie 
PID 


0.936 ±0.0071 




79,089 ± 379 


0.978 ±0.0024 


A 


1.62% ±0.007% 


FID 

/j.e 


0.928 ±0.0045 



Table 1 

Factors in the calculation of the B{K^l 
upper limit. 
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